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Abstract 

We study the fluctuation-induced, time-dependent force between two plates immersed in a fluid 
driven out of equilibrium mechanically by harmonic vibrations of one of the plates. Considering a 
simple Langevin dynamics for the fluid, we explicitly calculate the fluctuation-induced force acting 
on the plate at rest. The time-dependence of this force is characterized by a positive lag time 
with respect to the driving, indicating a finite speed of propagation of stress through the medium, 
reminiscent of waves. We obtain two distinctive contributions to the force, where one may be 
understood as directly emerging from the corresponding force in the static case, while the other is 
related to resonant dissipation in the cavity between the plates. 

PACS numbers: 64.70.qd,65.20.Dc,66.10.-x 
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I. INTRODUCTION 



A fundamental advance in the understanding of nature was the insight that physical forces 
between bodies, instead of operating at a distance, are generated by fields; the latter obeying 
their own dynamics, implying a finite speed of propagation of signals and causality 
Moreover, time-varying fields can sustain themselves in otherwise empty space to produce 
disembodied waves; exemplified by electromagnetic fields and waves, and gravitational fields. 
Gravitational waves are believed to be detected in the near future j^. 

Another force seemingly operating at a distance is the Casimir force. This force was first 
predicted by Casimir in 1948 for two parallel conducting plates in vacuum, separated by a 
distance L, for which he found an attractive force per unit area F/A = — 7r^/lc/(240L^) |3|. 
It can be understood as resulting from the modification of the quantum-mechanical zero- 
point fluctuations of the electromagnetic fields due to confining boundaries {4-^. In the 
last decade, high-precision measurements of the Casimir force have become available which 
confirm Casimir's prediction within a few per cent [7H9(; recent experiments demonstrate 
the possibility of using the Casimir force as an actuation force for movable elements in 



nanomechanical systems 



0,0. 



This development goes along with significant advances in 
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calculating the Casimir force for complex geometries and materials 

A force analogous to the electrodynamic Casimir force also occurs if the fluctuations 
of the confined medium are of thermal origin The thermal analog of the Casimir 

effect, referred to as critical Casimir effect, was first predicted by Fisher and de Gennes for 
the concentration fluctuations of a binary liquid mixture close to its critical demixing point 



confined by boundaries 
for this very system 



151 ]: recently, the critical Casimir effect was quantitatively confirmed 



16| . (For computational methods concerning the calculation of critical 



Casimir forces, see, e.g., references 
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The vast majority of work done on the Casimir effect, and fluctuation-induced forces 
in general, pertain to the equilibrium case. That is, the system is in its quantal ground 
state in case of the electrodynamic Casimir effect, or in thermodynamic equilibrium in case 
of the thermal analog. A number of recent experiments probe the Casimir force between 
moving components in nanomechanical systems j^, 10], and effects generated by moving 



boundaries have been studied, e.g., for Casimir force driven ratchets 



191 ]: however, the data 



are usually compared with predictions for the Casimir force obtained for systems at rest. 
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Figure 1. Two parallel plates separated by a varying distance L{t). Plate 1 is at rest while plate 2 
is vibrating parallel to the z-direction. The plates are immersed in a fluctuating medium described 
by the order parameter (j){r, t). The fluctuation-induced, time-dependent net force on plate 1 is the 
sum of the forces F^{t) and F- acting on opposite sides of the plate. 

corresponding to a quasi-static approximation. Distinct new effects occur if the fluctuating 
medium is driven out of equilibrium. In this case the observed effects become sensitive to 
the dynamics governing the fluctuation medium, which may lead to a better understanding 
of these systems, and may provide new control parameters to manipulate them 2[H24| . For 
example, the generalization of the electrodynamic Casimir effect to systems with moving 
boundaries, referred to as dynamic Casimir effect, exhibits friction of moving mirrors in 



vacuum and the creation of photons 



25N27l|. For the thermal analog, fluctuation-induced 



forces in non-equilibrium systems have been studied in the context of the Soret effect, which 
occurs in the presence of an external temperature gradient 211]. Fluctuation-induced forces 
lave also been obtained for macroscopic bodies immersed in mechanically driven systems 
2o| . granular fluids 22 1, and reaction-diffusion systems 23]. Recently it was shown that 
non-equilibrium fluctuations can induce self-forces on single, asymmetric objects, and may 
lead to a violation of the action- reaction principle between two objects 24]. 

In this work we consider a fluctuating medium driven out of equilibrium mechanically 
by a vibrating plate, and study the resulting time-dependent, fluctuation-induced force 
F(t) on a second plate at rest. We wish to elucidate the time-dependence of this force in 
view of the finite speed of propagation of signals in the fluctuating medium, and causality. 
Specifically, we consider two infinitely extended plates parallel to the xy-plane separated 
by a varying distance L{t) as shown in Fig.[Tl Plate 1 is at rest while plate 2 is vibrating 
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in ^-direction by some external driving. The plates are immersed in a medium undergoing 
thermal fluctuations with long-ranged correlations described by a scalar order parameter 
0(r, t). The order parameter is subject to Dirichlet boundary conditions = at the plates. 
As shown in Fig.[Tl F{t) is the sum of the forces F^{t) and F_ acting on opposite sides of 
the plate; F^{t) being the force acting on plate 1 from the side of the cavity, and F_ the 
(time-independent) force on the boundary surface of a semi-inflnite half-space flUed with the 
fluctuating medium. The net force F{t) = F^{t) + F^ is expected to be flnite and overall 
attractive, i.e., directed towards plate 2. 

Our presentation is organized as follows. In Sec. [IT] we introduce the Langevin dynamics 
of the order parameter 0(r, t) as a paradigmatic example for a non-equilibrium dynamics, 
and summarize our main results for the fluctuation-induced force F{t). In Sec. lIIII we discuss 
the calculation of F{t) using the stress tensor (Sec. lIII A|) and obtain F{t) to flrst order in the 



amplitude of the vibrations of plate 2 (Sec. lIII Bll . We find two distinct contributions to F(t), 
which can be attributed to real-valued poles (Sec. lIII Cl) and imaginary poles (Sec. lIII Dl) in 
the complex frequency plane, respectively, occurring in the calculation of F{t). We conclude 
in Sec.IVl 

II. MODEL AND MAIN RESULTS 

In the traditional case, both plates are at rest at a constant separation Lq (cf. Fig.[l]). 
The system is then in thermal equilibrium and the fluctuations of the order parameter are 
described by the statistical Boltzmann weight e~^^ with Gaussian Hamiltonian 

= l[d'r (V0)2 , (1) 



where (3 = l/{kBT) with the Boltzmann constant ks and the temperature T (assumed to 
)e const 



De constant). The fluctuation- induced force Fn on plate 1 per unit area A is found to be 



Fo _ J_{3)_kBT 
A~ Stt Lfj ' ^ ^ 

where the minus sign indicates that the force is attractive. Equation ([21) is a universal 
result, independent of the underlying dynamics of the fluctuating medium, as long as the 
thermodynamic equilibrium is described by Eq. ([1]). 
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We now turn to the case where plate 2 is vibrating parallel to the ^-direction, resulting 
in a time-dependent separation L{t) to plate 1. The time-dependent boundary conditions 
for the order parameter 0(r, t) then drive the system out of equilibrium. Locally, the order 
parameter will relax back to equilibrium according to the dynamics of the medium; in this 
work, we consider an overdamped dynamics described by the Langevin equation 

7^0(r,t) = VV(r,t)+r/(r,t) (3) 

where 7 is the friction coefficient. The random force ri{r,t) is assumed to have zero mean 
and to obey the fluctuation-dissipation relation 

{7^{r,t)v{r',t')) = 2^kBT5^'\v - v')5{t - t') (4) 

where the brackets ( ) denote a local, stochastic average and 5'^^^ is the delta function in 3 
dimensions. 

Figures [2] - m summarize our main results for the case that the external driving L{t) is 
given by harmonic oscillations 

L{t) = Lo + acos(u;ot) (5) 

with amplitude a and frequency ujq. Our results for F{t) hold to first order in a and can be 
cast in the form 



-^0 



0{a') (6) 



where Fq from Eq. ([2]) is the force for a constant separation Lq. The dimensionless parameter 

n = uoiLl (7) 

characterizes the strength of the friction coefficient 7 of the fluctuating medium (cf. Eq. ([3])). 
Equation ([6]) implies that the dimensionless function f{t,Q) is normalized such that / = 1 
for uq = 0, i.e., = 0. It can be represented as 

/(t, n) = A cos(cjot = A cos[u;o(t - to)] (8) 

in terms of an amplitude A{Q) and a phase shift ip{Q). The second equation in Eq. 
holds if ip is proportional to the driving frequency uq, i.e., (p{0,) = uJoto{Q), where the lag 
time to corresponds to the time delay between the source (driving L{t) of plate 2) and the 
resulting response (force F(t) at plate 1). For the relation (f{^l) = co'oto(^) used to obtain 
the second equation in Eq. ([H]) it is understood that to(^) depends on Q only weakly so that, 
in particular, to{0) is flnite. 
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Figure 2. Amplitude A of f{t,^l) as a function of (see Eq. ([8])). Shown are results for F{t) (black 
line) and for the contributions to F{t) according to Eq. ((33]) (red line) and Eq. (f35|l (green line), 
respectively. 




Figure 3. (a) Lag time to in terms of the dimensionless combination r = to/(7-f'o) ^ function of 
O (see Eq. ([S])). Shown are results for F{t) (black line) and for the contribution to F{t) according 
to Eq. ([33|) (red line). In both cases, the dependence t{Q), i.e., to{u!o), is fairly weak, so that the 
interpretation of to as a lag time is justified, (b) Phase shift (/9 as a function of (see Eq. ([8])) for 
the contribution to F(t) according to Eq. ([35|) . In this case, to = ^/^o strongly depends on loq. 

Figure [2] shows the amplitude A{fl) of the function /(t, Q) according to Eq. ([8]) (black 
line). The contribution to F{t) according to Eqs. (1331) and fl35|l . corresponding to real- valued 
poles and imaginary poles in the complex frequency plane occurring in the calculation of 
F{t), are also shown (red and green lines, respectively); cf. Secs. lIII Cl and HTl DI below. Figure 

shows the lag time to(fi) according to Eq. ([H]) in terms of the dimensionless combination 
r = to/(7-^o)- Shown are results for F(t) (black line) and for the contribution to F(t) 
according to Eq. fl33|l (red line). In both cases, the dependence t{Q), i.e., tol^^o), is fairly 
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Figure 4. c{t) = F{t)/FQ according to Eq. ([6]) for cjq = 27rs~^ and Q = W (black line). The function 
cos(woO) corresponding to the oscillating part of L{t) with a = 1 (cf. Eq. ([5])) is also shown (blue 
line). The lag time for these parameters is to = t{Q,)Q/ujo ~ 0.15s. 

weak, so that the interpretation of to as a lag time is justified. Thus, the form of f{t, Q) in 
Eq. (IH]), in conjunction with an approximately constant lag time to, indicates that the stress 
generated locally at the vibrating plate 2 is carried through the medium, according to its 
diffusive dynamics, until it arrives at plate 1 after a time to ~ 7-^o- This indicates that the 
fluctuation-induced force on plate 1 is indeed not operating at a distance, but generated by 
a local field, which for the present system is presumably related to the local stress in the 
fluctuating medium between the plates. In contrast, Fig.[3]D shows that for the contribution 
to F{t) according to Eq. ( l35l) . the phase shift (f rather than to is fairly constant, which implies 
that to = ^p/ujq strongly depends on Uo- Note, however, that for small cjq, i.e., small Q, the 
corresponding contribution to F{t) is suppressed by a vanishing amplitude A{Q) (green line 
in Fig. [2]). For illustration. Fig. H] shows the time-dependent force F{t), normalized by its 
value Fq for uq = 0, for the values Uo = 27rs~^ and 1] = 10 (black line). The function 
cos(ci;o^)) corresponding to the oscillating part of L{t) with a = 1 (cf. Eq. ([5])), is also shown 
(blue line). The lag time for these parameters is to = t{Q)Q/ujo — 0.15s. 
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III. METHOD 



A. Calculation of F{t) using the stress tensor 

The force per unit area acting on plate 1 from the side of the cavity can be expressed as 
F+{t)/A = \im{Tzz{r\\, z,t)) where ry = {x,y) are the components of r parallel to the plate 
and = \ idz(f))'^ — | [{dx(f))'^ + (dy(t>)^ ] is the 2;^-component of the stress tensor, assuming 
the system is locally in equilibrium 



281]. Similarly, the force per unit area acting on the 
other side of plate 1 obtains as F^/A = — lim(T22(r|| , z, t))L=oo where Tzz is again evaluated 

2— >0 

in the cavity between the plates but for the limit L — > oo (cf. Fig.[l]). The net force per unit 
area on plate 1 obtains as 

Fit) 

—— = lini(T,^(r||,2;,t)) - lim(T^^(r||, z, t))i=oo (9) 

/i 2-^0 2— +0 

where, using the Dirichlet boundary condition = at the plates, 

lim(T,,(r||,z,t)) = l lim dzdz'{(p{T\\, z,t)^iT\\, z' ,t)) . (10) 

2^0 Z 2,2'^0 



To calculate the two-point correlation function of on the right-hand side of Eq. (11 01) we 
note that the solution 0(r, t) of Eq. ([3]) can be expressed as 

oo 

^{v,t)= j dt' j £r'G{r,t;r',t')r]{v',t') (11) 

-oo V{t') 

where V{t') = A ■ L{t') is the volume of the cavity at time t' and the Green's function 
G(r, t; r', t') is defined as the solution of 



- V,^) G(r, t- r', t') = 5(3) (r - r')5(t - t') 



(12) 



subject to the boundary condition G(r, t; r', t') = whenever r or r' is located on the surface 
of one of the plates at time t or t', respectively. In addition, G{r,t] r',t') = for t' > t hj 
causality. Thus, 0(r, t) can be expressed as a linear superposition of contributions from 
the source ri{r',t') at times t' < t and positions r' G V{t'), carried forward in time by the 
propagator G(r, t; r', t'). Using Eqs. (fTTl) and (jl]), the two-point correlation function of 
obtains as 

00 

(0(r,t)0(r',t')) = 27A;bT j ds j d^x G{r,t;x, s)G{r' ,t';x, s) . (13) 

-00 V{s) 
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In the present set-up, the system is translationally invariant in xy-direction at any time t, 
whereas translation invariance in time is broken due to the varying separation L{t) between 
the plates. Thus, introducing the partial Fourier transform of G as 



G(r||, z, t; rjj, z', t' 



d'^P ip-(r||-r;) f duj 



(2vr) 



-iu){t-t') 



2-K 



g{z,z';uj,p, t' 



(14) 



the function g depends explicitly on one of the time coordinates in G, say, t'. Using Eqs. (1131) . 



M|) we find for the expression in Eq. ([TO 



3Q| 



hm{T^^{r\\,z,t)) = 'jhsT 



(Pp f duj 



(27r)2 J 2tt 

-oo 



dCu{C,uj,p,t)u*{C,uj,p,t) (15) 



where 



d 

u{C,uj,p,t) = —g{z,C,^^,P,t) 



z=0 



(16) 



For given propagator G, hence function u, F{t)/A can be calculated using Eqs. ([9]) and (fT5l) . 



B. Calculation of the propagator G 

The remaining task is to calculate the propagator G(r, t; r', t') solving Eq. (fT2|) subject 
to the time- dependent boundary conditions due to the vibrating plate 2. This problem 
can be solved, for general modulations of the plate(s) in space and time, by the method 



developed in reference 311] • For the present set-up, we find for the partial Fourier transform 



Q{z,t] z',t']p) of G (i.e., transforming the spatial coordinates r||, rji parallel to the plates as 



in Eq. ( fT4l) but keeping 



for ease of notation) 3l| 



the time coordinates t, t'; in what follows, we omit the argument p 



g{z,t;z',t') =g{z,t;z',t') - J dr J da (17) 

— oo — oo 

g[z, t; L{t), t]M{t, a) g[L{a), a; z\ t'] 
where Q is the propagator in the half-space 2; > bounded by a Dirichlet surface at 2; = 



32l |. The kernel M. is defined by 

00 

j daM{r,a)g[L{a),a;L{t),t\=5{T-t) . (18) 



In this work, we consider small variations of the separation between the plates about a 
mean separation Lq, i.e., 

L{t) =Lo + h{t) . (19) 
Our results hold to first order in h. To this end, we insert Eq. ( fT9l) in Eq. ( fT7|) and expand 



everything to first order in h 33|. This results in expansions g = Qo + Qi + 0{h'^) and 
u = Uq + Ui + 0{h'^) of the functions g and u from Eqs. fll4|] . fITB]) in powers of h. Equations 
(E]), f|T5l) then yield the corresponding contributions to F{t)/A. 

Let us first consider the leading order, i.e., /i = and L{t) = Lq. Using Eq. (fT7I) and 
transforming to tu-space as in Eq. f ll4p we find (omitting the arguments p and u for ease of 
notation) 

go{z, z') = giz, z') - g{z, Lq) Mq ^(Lq, z') (20) 

where g{z,z') = |^e~*5l^~^'l — e"'^'^^"*"^'^] /{2Q) with Q = a/p^ — i7Co' from Eq. (IMl) below and 
Mo = [^(Lo,Lo)]-^ = 2g/[l - exp(-2gLo)]. Thus, 

gsmh(gLo) 



and, using Eq. ([16 



(92;'^ 



sinh[Q(Lo - 0] 



UoiC) = Trgo{z, C) = .^'"r ^ ■ (22) 



2=0 sinh (QI. 



Using Eqs. (jH]), (H^H) . (12^ we thus obtain to leading order [21 1 



4) J 



oo 

[(%[ [coth(QLo) - 1] - P [coth(PLo) - 1]) (23a) 



Fo _ ksT f d^p f du) 1 
~A ~ 2 



-oo 



knT f d^P r , / N n / , V 

' ^ p [coth(pLo) - 1] . (23b) 



2 J (27r)2 

The integral in Eq. (I23bl) is finite and yields Eq. ([2]). In Eq. fl23al) we use 

Q{uj,p) = ^/j^^^iyu , P{uj,p) = ^/p^ + i-fu , (24) 

so that P = Q* if u is real. Integrations over uj as in Eq. fl23ap are readily computed by 
contour integration in the complex cu-plane. In Eq. f l23al) and throughout this work we use the 
convention that in cj-integrations we integrate above the pole in u; this can be accomplished 
by the replacement u ^ u + ie in the denominator of the integrand in Eq. fl23al) . The 
limit 5 — i> in final results is always understood. Note that this prescription introduces 
a positive time direction and ensures causality. Q{uj) has a branch cut along the negative 
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imaginary axis 70; = —i{p'^ + r), r > 0, whereas -P(ci;) has a branch cut along the positive 
imaginary axis 7C1J = z(p^ + r), r > 0. The integral over uj in Eq. fl23al) has two contributions. 
For the contribution involving Q [coth((5i^o) ~ l]? the contour integral can be closed in the 
upper complex cj-plane (thus avoiding the branch cut of Q), where this term has no poles, 
so that the contribution from this term vanishes. Likewise, for the contribution involving 
P [coth(PLo) — 1], the contour integral can be closed in the lower complex cj-plane (avoiding 
the branch cut of P), where, in turn, this term has no poles. The single pole sX uj = —is 
in the lower complex cu-plane then yields the expression in Eq. (123bl) : cp. Fig. [5^ in Sec. lIII Cl 
with ujq = 0. 

We now turn to the contribution to F{t)/A to first order in h. Using the expansion 

uiC;uj,p,t) = Mo(C;t^,p) + Mi(C;t^,p,i) + O {h"^) (25) 

in Eq. (fTSjl . with u from Eq. (fT6|) and Uq from Eq. (1221) . we find for general h{t) 

oo 

F.jt) ^ kBT f d'p 
A 

where 



/(^/ ^[/(^'P'^) + r(^'P'^)] (26) 



-oo 



Q h 1 
f{uj,P,t) = . o 



Q 



smh{QLo) iuj ysmhi^QLo) sinh(PLo) 



(27) 



h 



The symbol o denotes a convolution of two functions a{uo)^ b{uj) involving an insertion of 



oo 

(a o b){u, t') = d{u) I — e-^(— ')*'/i(cu - uo')h{uj') . (28) 

27r 

For functions a{t,t'), b{t,t'), the expression (a o b){uj,t') is the representation in u;-space of 
c(t,t') := JZ,dsa{t,s)his)bis,t'y, i.e., c{t,t') = J^^^ exp[-tujit - t')]id o b){uj,t').The 



functions d{uj), b{uj) are the representations in c<j-space of a(t, s), b{s,t'), respectively |34l |. 

For the special case that plate 2 is vibrating with harmonic oscillations of amplitude a 
and frequency uq (cf. Eqs. (jS]), ( fT9l) ). i.e., 

h{t) = acos{LUot) , (29) 

we obtain h^u) = an [6{u - ujq) + 6{uj + ujq)]. The integral ^{f + f*) in Eq. ([26]) decays 
into two contributions corresponding to the terms in square brackets on the right-hand side 
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of Eq.([27D: 

OO 

/duj 
-^[f{u,p,t) + r{uj,p,t)] = ^QQM + ^Qp{p,t) , (30) 

— oo 

where the subscripts QQ and QP indicate the contributions from the first and second term 
in square brackets of Eq. (1271) . respectively. In what follows we show that these two terms 
yield distinct contributions to F{t) corresponding to real-valued and imaginary poles in the 
complex cu-plane. 



C. Real- valued frequency poles: lag time to 



For the first contribution in Eq. fl30l) we find [30|, l35 | 

oo 

u{uj)u{(jJ — OJq) v{ijJ)v{ijJ + OJq) 

oj — ojq + ie uj + ujq + ie 



+ c.c. (31) 



where 



smh((yLo) smh(PLo) 
with Q, P from Eq. (1241) . Computing the right-hand side of Eq. ( I3T1) by contour integration 

in the complex tu-plane, the contributions from the two terms in square brackets in the 
integrand are analyzed along similar lines as discussed below Eq. (124|) . Thus, for the first 
term in square brackets, the contour integral can be closed in the upper complex cj-plane, 
where u{uo) has no poles, so that the contribution from this term vanishes. For the second 
term in square brackets, the contour integral can be closed in the lower complex cj-plane, 
where v{uj) has no poles. The only contribution from this term is from the single pole at 
uj = —uq — ie; see Fig.|5^. Thus, including the contribution from the complex conjugate in 
Eq. f l3T]) . we obtain 

The corresponding contribution to Fi{t)/A is given by J -^^(^QQ{p,t) (see Eqs. ( 126|1 
and (l30l)). In the static case, where cuq = and h{t) = a in Eq. ([5]), this result can also 
be obtained directly from Eq. (I23bl) by replacing Lq with Lq + a and expanding to first 
order in a. For finite ujq, Eq. (133|) emerges from the static case by a shift from cuq = to 
a finite value of ujq. This shift may be understood in terms of a transition from stationary 
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modes (standing waves) in the cavity in the static case to modes with a time-dependence 
~ exp(zci;ot), reminiscent of travehng waves, in response to the oscillating plate 2. The time- 
dependence of these modes carries over to the fluctuation-induced force on plate 1 (cf. Fig.[T]). 
The picture of traveling waves in the fluctuating medium with a finite speed of propagation 
(diffusion) is consistent with the presence of a lag time to ~ 7-^o ^'i- dS]); where to > by 
causality; cf. Fig. [3^1 and the related discussion in Sec.HTl 



D. Imaginary frequency poles: Resonant dissipation 



For the second contribution in Eq. (I30l) we find |35| 



oo 



{u — ujq + ie){ijj + is) 



+ c.c. (34) 



The contour integral over uj can be closed either in the upper or the lower complex u;-plane, 
yielding identical results; the contributions from the poles aX uj = uoq — ie and uo = —ie 
in the lower complex cu-plane cancel. Closing the contour integral in the lower plane, the 
integral picks up contributions from the imaginary poles 'jUn = —i{p^ + k^) of u{uj), where 
kn = niT / Lq and n > 1 is a positive integer. Note that u{uj) has a branch cut along 
the negative imaginary axis on which the poles uOn are located (cf. the related discussion 
below Eq. fl24|) ): however, this branch cut may be cured using the identity Q/ smh{QLQ) = 



R/ sm{RLo), with R{u,p) = ^/i'yuj — p^, which holds close to the negative imaginary axis. 
The expression R/ sin(i?Lo) is analytic in the lower complex cu-plane with isolated poles at 
ujn] see Fig.Eb. Summing over the residues of these poles yields 



Lq ^ (r + ^„)(r + /Sn-^7^o) 

The corresponding contribution to Fi{t)/A is given by / ■^^^$Qp(p, t). Note that 
^Qp{p,t) is proportional to ujq, which implies that this term is absent in the static case 
uq = and solely generated by the fact that the system is driven out of equilibrium by the 
oscillating plate. The imaginary poles 70;^ = + k"^) leading to Eq. fl35|) are related to 

the dissipation in the medium. The poles Un thus correspond to resonant dissipation, where 
the spectrum of resonance frequencies -|- n'^ir'^ /{'jL'^) is continuous due to the presence 
of the continuous in-plane wave number p 36 1. 
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a) 



b) 



complex 
CO — plane 
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Figure 5. (a) Contour integration in the complex cj-plane for the second term in square brackets 
in Eq. (j3ip . The only contribution from this term is from the pole at a; = —ujq — ie indicated by 
the red star. The blue lines indicate branch cuts of v{uj)v{u> + ujq). (b) Contour integration in 
the complex w-plane in Eq. (j34p . The blue lines indicate branch cuts of of u{uj)v{u> — ujq). The 
contributions from the poles at to = ujq — ie and uj = —ie cancel (see text). 

IV. CONCLUDING REMARKS 

We have studied the time-dependent, fluctuation-induced force F{t) on a plate at rest 
generated by a second plate with harmonic oscillations (cf. FigH]). Our main results, valid to 
flrst order in the amplitude a of the oscillations (cf. Eq. (^) and summarized in Figs.[2]-|31 
indicate that the fluctuation-induced force is carried through the medium from one plate to 
the other with a flnite speed of propagation (diffusion). We flnd two distinct contributions 
to F{t), related to real- valued poles and imaginary poles in the complex frequency plane, 
resulting in a flnite lag time to in Eq. ([8]) and resonant dissipation (Secs. HTTCl and Ull Dl) . In 
this work we consider a scalar order parameter 0(r, t) with overdamped dynamics described 
by the Langevin equation (cf. Eq. ([3])). However, our approach can be readily applied to 
other dynamical systems. In particular, it would be interesting to extend this study to the 
time dependence of the electro dynamic Casimir force. 
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